class Riemann_Kernels
{
public:
  inline Riemann_Kernels(const field<mat> in_riemann_points,  int in_dim);
  
  inline mat Stein( std::string tosave_in);
  inline mat Stein_faster( std::string tosave_in);

  
  
private:
  const field<mat> riemann_points;
  int N_points;
  int dim;
 
  
};